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Abstract 


The inelastic hadronic cross section in proton-lead collisions at a centre-of-mass en¬ 
ergy per nucleon pair of 5.02 TeV is measured with the CMS detector at the LHC. 
The data sample, corresponding to an integrated luminosity of £ = 12.6 ± 0.4nb^^, 
has been collected with an unbiased trigger for inclusive particle production. The 
cross section is obtained from the measured number of proton-lead collisions with 
hadronic activity produced in the pseudorapidity ranges 3 < rj < 5 and/or 
—5<rj< —3, corrected for photon-induced contributions, experimental accep¬ 
tance, and other instrumental effects. The inelastic cross section is measured to be 
Hnei(pPt>) = 2061 ± 3 (stat) ± 34 (syst) ± 72 (lumi) mb. Various Monte Carlo genera¬ 
tors, commonly used in heavy ion and cosmic ray physics, are found to reproduce the 
data within uncertainties. The value of (7jj.iei(pPb) is compatible with that expected 
from the proton-proton cross section at 5.02 TeV scaled up within a simple Glauber 
approach to account for multiple scatterings in the lead nucleus, indicating that fur¬ 
ther net nuclear corrections are small. 


Published in Physics Letters B as doi : 10.1016/j.physletb. 2016. 06. 02 7. 


© 2016 CERN for the benefit of the CMS Collaboration. CC-BY-3.0 license 


*See Appendix|^for the list of collaboration members 










1 


1 Introduction 

The measurement of the inelastic cross section in proton-lead collisions, c7'inei(pPb), at a centre- 
of-mass energy per nucleon pair of 5.02 TeV performed by the CMS experiment at the CERN 
LHC is presented. The inelastic cross section (also called "particle-production" [IJ or "absorp¬ 
tion" |j2l cross section in previous studies) is defined to include all hadronic events, including 
contributions from diffractive processes, except those from the quasi-elastic excitation of the 
lead nucleus—estimated to amount to about 100 mb for the pPb system |@. Inelastic electro¬ 
magnetic (photon-proton) collisions are also excluded from the measurement. 

While being one of the most inclusive observables in hadronic collisions, the inelastic cross sec¬ 
tion is one of the least theoretically accessible quantities, as it cannot be determined from first- 
principles calculations of the theory of the strong interaction, quantum chromodynamics. In 
proton-proton (pp) and nucleus-nucleus collisions at the LHC, particles produced in hadronic 
interactions come mostly from the hadronisation of quarks and gluons, either produced in 
semi-hard scatterings ("minijets") |4] or emitted at very forward rapidities from "spectator" 
partons, as well as from soft diffractive processes in "peripheral" interactions. Lrom the mea¬ 
sured inelastic proton-proton (or nucleon-nucleon) cross section at a given collision energy, one 
can theoretically derive the corresponding proton-nucleus and nucleus-nucleus cross sections 
by means of Glauber JSl |6l or Gribov-Regge IZJ multiple-scattering approaches that take into 
account the known transverse matter profile of nuclei. Key quantities for the experimental com¬ 
parison between nucleus-nucleus and pp collisions—such as the nuclear overlap function, the 
number of nucleon-nucleon collisions and of participant nucleons [Sj 0 —are also commonly 
computed through such approaches. Validating the Glauber and Gribov-Regge predictions 
with proton-nucleus collisions at LHG energies has important implications beyond collider 
physics. Such approaches constitute crucial ingredients in the Monte Garlo modelling of cos¬ 
mic ray air showers at the highest energies IfTOl , for which the inelastic cross sections measured 
in the laboratory must be extrapolated over a wide energy range. In fact, the inelastic proton-air 
(mostly proton-nitrogen and proton-oxygen) cross section introduces one of the largest uncer¬ 
tainties for air shower simulations iiiiiia. 

The Glauber multiple-collision model, based on the eikonal limit (i.e. straight-line trajectories 
of the colliding nucleons), is the simplest and most economical approach often used to derive 
inclusive proton-nucleus quantities from the pp cross sections and, vice versa, to obtain pp 
cross sections from the cosmic ray measurements IIT3ll . However, some of the approximations 
applied in the model—foremost the absence of short-range nucleon correlations I1T4I and of 
inelastic screening flSl —impact the computed cross section values. This is observed for fixed- 
target proton-carbon data IIT6] - (20l and estimated for collider ||2Tl|22l as well as ultra-high cosmic 
ray IT^ energies, where corrections to the proton-air cross section of the order of 10% have been 
obtained. Short-range correlations increase the number of nucleon-nucleon collisions at small 
impact parameters yielding a larger nucleus-nucleus cross section. On the other hand, screen¬ 
ing affects the number of nucleons that are diffractively excited in the multiple collisions but 
revert back to their ground state before the scattering process is completed, thereby reducing 
the nuclear cross section. Different implementations of such effects exist in the current hadronic 
interaction models llT5ll23ti29| . A measurement of cT-nei in pPb collisions at the LHG can test if 
the precision of the standard Glauber calculation is sufficient, and at which energies corrections 
to the Glauber approach may become relevant. 
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3 Event selection and analysis 


2 Experimental setup and Monte Carlo simulations 

The measurement presented here is based on pPb data taken with the CMS experiment at the 
LHC at the beginning of 2013. A detailed description of the apparatus can be found in If3^ i|. The 
main detector used in this analysis is the hadron forward (HP) calorimeter that covers the pseu¬ 
dorapidity interval 3 < |7/| < 5. The calorimeter is composed of quartz fibres in a steel matrix 
with a 0.175x0.175 segmentation in the azimuthal angle (p (in radians) and pseudorapidity rj. 
The quartz fibres pick up the Cherenkov light produced by the charged component of showers. 
This light is then measured by photodetector tubes. The hadronic and electromagnetic signals 
of each segment, as derived from fibres of two different lengths and depths, are combined to 
form a tower signal. 

The data used in this analysis comprise an integrated luminosity of C = 12. 6 ± 0.4 nb This 
dataset combines the integrated luminosities of the two possible directions of the proton and 
lead beams: 5.0 ± 0.2 nb^^ and 7.6 ± 0.3 nb~^, for the proton beam going respectively in the 
clockwise (negative t]) and anticlockwise (positive ?/) direction. The events are collected using 
an unbiased trigger, only requiring the presence of both beams in the interaction point, as 
determined by the "Beam Pickup Timing for the experiments" (BPTX) devices. Detector noise 
is studied with events that are randomly read out in the absence of both beams in the detector. 
The luminosity determination technique was calibrated by means of a van der Meer scan iMlI 
for both beam directions independently, with an uncertainty of 3.5% 1321 . 

A Monte Carlo event simulation based on a Ge ANT4 detector description Il33l is used to model 
the experimental response and derive the reconstruction efficiencies. Different event genera¬ 
tors are used to simulate hadronic proton-nucleus collisions. Three models are based on the 
Gribov-Regge formalism: DPMJET 3.06 IIMI . EPOS-LHC Il25l , and QGSJETll-04 Il26l : and a fourth 
one is based on a minijet+Glauber approach: HIJING 1.383 11351 . In addition, particle produc¬ 
tion from photon-proton (yp) interactions in "ultraperipheral" collisions, at impact parame¬ 
ters larger than the sum of proton and lead radii, needs to be taken into account ||3^ . Given 
the large Pb ion charge, and the associated large "equivalent photon flux" of its electromag¬ 
netic field O/ inelastic photon-proton collisions result in a non-negligible particle production 
contribution. Pure photon-photon interactions, mostly producing exclusive electron-positron 
pairs, and photon-nucleus interactions (where the photon emitted from the proton collides 
with the Pb ion) have orders-of-magnitude smaller visible cross sections and are neglected. 
Photon-proton processes are generated with the starlight programme 1371 combined either 
with DPMJET 3.05 or PYTHIA 6.4.26 1381 . 

3 Event selection and analysis 

In this analysis three types of cross sections are measured: (i) Gobs after removal of noise and cor¬ 
rection for pileup, (ii) Gvis after further removal of electromagnetic contributions and translation 
into a hadron-level quantity, and (iii) Gjoei including the final extrapolation to the total inelastic 
hadronic cross section. Two different approaches are used to determine the number of inelas¬ 
tic events: (1) a single-arm event selection that requires a localised calorimetric energy signal 
above a given threshold in the HP detector either at positive or negative pseudorapidities, and 
(2) a double-arm event selection that requires a localised signal above threshold in both HP 
detectors. The advantage of using these two event selections is that they have very different 
sensitivities to diffractive and photon-proton events as well as to detector noise. Denoting by 
£hf+ (Trf-) the highest energy measured in an HP tower at positive (negative) pseudorapidity. 



3 


an event is tagged as a candidate for an inelastic collision if it has a value of 

max ( Ehf+/ Erf- ) for single-arm selection 


Erf — 


min ( Ehf+/ Ehf- ) for double-arm selection 


( 1 ) 


above a given threshold. 
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Figure 1: Distribution of the energy deposited in the HF calorimeter (Ehf) for the single-arm 
(left) and double-arm (right) event selections. The data sample, shown exemplarily for one 
period with stable run conditions, comprises 1.31 nb^^ recorded with an unbiased trigger. The 
contribution from noise is obtained from a random trigger normalised to the same number of 
triggers as that in the collision data. The average number of photon-proton processes simulated 
with STARLIGHT+DPMJET and STARLIGHT+PYTHIA is treated as background and stacked on 
top. Four hadronic interaction models (EPOS, DPMJET, HIJING, and QGSJETIl) are overlaid and 
normalised to the number of data events with Ehf > 10 GeV, where the contribution from the 
background is small. The vertical line represents the threshold energy of 8 GeV (4 GeV) for the 
single-arm (double-arm) selection used in this analysis. 


The observed distribution of Ehf is well reproduced by the combined hadronic inelastic, photon- 
proton, and detector noise contributions as shown in the left (right) panel of Fig.j^for the single¬ 
arm (double-arm) selection. The size of the various contributions to the HF energy deposition 
is determined from data and simulations. The signal is identified as that coming from hadronic 
collisions whereas the backgrounds arise from electromagnetic photon-proton interactions and 
detector noise. The expected number of photon-proton collisions is Nyp = f^pCr^pC, where /yp 
is the fraction of simulated photon-proton events passing the selection and CFyp is the predicted 
STARLIGHT cross section. The number of misidentified events produced by electronic noise in 
the detector is Nnoise = Nfnoiser where /noise is the fraction of events read out randomly in the 
absence of beams that pass the selection criteria, and N is the number of events recorded with 
the unbiased trigger. The estimate of Nnoise includes Nobs+noise = Nobs/noise events that contain 
also an observed inelastic collision, where Nobs is the number of observed inelastic events. The 
double-counted events are explicitly subtracted from Nnoise- The uncertainty on Nnoise is de¬ 
rived from variations in different data-takmg periods. The background induced by beam-gas 
collisions is found to be negligible deduced from the fraction of events selected with the trigger 
indicating the presence of a single beam in the interaction point. 
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3 Event selection and analysis 


Of the number of inelastic hadronic collisions, NmeU the ones that are observed by the de¬ 
tector and pass the event selection are defined as Nhad- The purity of the event selection is 
Nhad/ (Nhad + N^p + Nnoise)/ and the acceptance is given by the ratio Cacc = Nhad/Ninei- Both 
the purity and the acceptance depend on the energy threshold used for the selection. Higher 
purity is achieved for the double-arm selection, since photon-proton interactions lead to a typ¬ 
ical final state where most of the secondary products are asymmetrically emitted towards the 
direction of the proton beam. Noise events are also suppressed by the coincidence requirement. 
The acceptance is in general smaller for the double-arm selection due to the smaller chance of 
selecting diffractive events characterised by large rapidity gaps devoid of activity in one or 
both HP sides. 

The dependence of Cacc on the HP tower energy threshold is shown in Pig. Por the single¬ 
arm selection the working point is chosen to be Phf > 8 GeV. This value is the result of a 
compromise between acceptance (about 93-94%) and contamination, while the probability to 
have a tower above the threshold does not depend much on the beam direction. The double¬ 
arm selection uses £hf > 4 GeV yielding 99% purity and 91% acceptance. The value e^cc for 
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Pigure 2: Acceptance versus purity of the two event selections, as derived from the EPOS and 
QGSJETll generators. The symbols indicate different values of the Phf thresholds. The chosen 
thresholds are marked with squares. 

a specific £hf threshold is determined by averaging over the results of the EPOS and QGSJETll 
models. The results of HIJING and dpmjet, which do not include nuclear effects for diffraction, 
are not considered for this purpose. Indeed, we have verified that both latter models are unable 
to describe the very-forward energy spectra measured with the GASTOR detector (—6.6 < 
rj < —5.2) in events with large rapidity gaps, which are particularly sensitive to diffractive 
interactions. 

The uncertainties on the Cacc and Nyp values are estimated from the maximum absolute dif¬ 
ferences obtained from the results of different event generators, averaged over a wide £hf 
interval between 2 and 10 GeV. The uncertainties on Cacc are 0.005(0.014) and of N^p/ C are 
11 mb (0.05 mb) for the single-arm (double-arm) event selections. 
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In this analysis no vertex reconstruction is performed and the impact of contributions from ad¬ 
ditional pileup (PU) collisions recorded in any given event is consistently evaluated with the 
HP detector. The number of simultaneous collisions is Poisson-distributed with an expectation 
value corresponding to the interaction probability A. If one collision is selected with probability 
6acc/ then i simultaneous collisions are selected with probability Pi « 1 — (1 — Cacc)*- The ap¬ 
proximation assumed in the equation, which does not account for energy deposits of multiple 

events in the HP towers, was verified to be valid by means of a toy Monte Carlo simulation. 

00 

The number of collisions is then corrected with the factor /pu = ^accT/ YL Pi Poisson(z; A). The 

i=l 

interaction probability A, which amounts to 2-8% depending on the data-taking period, is cal¬ 
culated recursively from the ratio of the number of inelastic events to the number of unbiased 
triggers. The pileup correction increases the measured cross section by 2% for both event selec¬ 
tions, and introduces an uncertainty on the final pPb cross section that is smaller than 0.1%. 

To facilitate the direct comparison of the results to model predictions, detector level quanti¬ 
ties, such as £hf/ are translated to hadron-level quantities. For this purpose, Phf is defined 
equivalently to Eq. Q but replacing Ehf by the largest absolute value among the momenta, \p\, 
of all generated final-state particles (with lifetimes above 1 cm/c), within the pseudorapidity 
intervals of the HP calorimeters (3 < | //1 < 5), excluding muons and neutrinos. A correction 
factor Cvis/ obtained from simulations, is used to translate the measured cross section into a 
hadron-level quantity, defined by the ratio of the number of visible events, which fulfil a given 
requirement on phF/ to the number of observed events, which pass the selection on Ehf- Thus, 
Cvis is larger than unity for requiring Phf > 0, but will approach zero for very high thresholds. 
The threshold can be chosen freely, and for the present analysis the requirement on the minimal 
value of Phf is chosen such that the fractions of events passing this selection and passing that 
on Ehf are equal. The factor Cvis then becomes equal to unity and has no numerical effect on the 
central value of the derived cross section. This procedure leads to the choice of selecting events 
that fulfil the requirement Phf > 21.3GeV/c (11.3GeV/c) for the single-arm (double-arm) anal¬ 
ysis. For the chosen thresholds, the mean of the Cvis values of all four hadronic interaction 
models is unity and the slight dependence on models is taken into account as a systematic 
uncertainty on Cvis equal to the standard deviation of the four values. 

The values of the acceptance, backgrounds, and correction factors are summarised in Table 

Table 1: Gentral values and uncertainties for the two event selections for noise cross section 
contribution (Nnoise/ P) arid the fraction of noise events (/noise) as derived from data. Addition¬ 
ally, the quantities acceptance (Cacc)/ electromagnetic cross section contribution (Njp/C), and 
hadron-level correction factor (Cvis) as derived from simulations are listed. 


Selection 

Nnoise /P [Pob] 

/noise 

^^acc 

Nypf C [mb] 

Cvis 

Single-arm 

102 ± 25 

(2.0 ±0.5) X 10-^ 

0.939 ± 0.005 

63 ±11 

1.000 ±0.004 

Double-arm 

9±3 

(1.8 ±0.8) X 10-4 

0.910 ± 0.014 

0.33 ± 0.05 

1.000 ±0.002 


The number of observed inelastic events. Nobs/ is derived from the number of events passing the 
event selection, Ngei, and is corrected for noise (Nnoise)/ double counting (Nobs+noise)/ and pileup 
(/pu). Dividing this number by the integrated luminosity yields the observed cross section: 


G)bs 


Nobs 

~C~ 


(Nsel - Nnoise 3“ Nobs+noise) 


/pu 

C ■ 


( 2 ) 


Using the relation Nobs+noise = Nobs/noise one obtains 

_ 1 Nggi Nnoise 

E1//PU-/noise 


( 3 ) 
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3 Event selection and analysis 


The visible cross section for hadronic collisions is derived by subtracting the photon-proton 
contamination and applying the correction factor Cvis- Its numerical value is, by definition, 
equal to the part of the observed cross section related to hadronic collisions: 


_ 1 hfsel hinoise N7P ^ 

^ 1 / /pu /noise 


(4) 


The inelastic cross section is obtained by correcting for the limited detector acceptance (Cacc): 


^inel 


1 hsfnoise 1 

^ 1 / /pU /noise ^acc 


(5) 


The ratio of the visible hadronic cross section obtained with the single-arm selection to the one 
obtained with the double-arm selection is sensitive to the fraction of diffractive pPb events. 
The measured value of this ratio allows one to constrain the diffractive cross section, in 
the models. In order to be compatible within 2 standard deviations of the data, the EPOS diffrac¬ 
tive cross section cannot be scaled up or down by more than ±13% from its default value, while 
for QGSJETll those limits are ±20%. This propagates into an eacc(Pdiff) uncertainty on Pinei, con¬ 
servatively assumed to be symmetric, of 0.8% (1.1%). For this and the following uncertainties, 
the first number is related to the single-arm selection and the bracketed one to the double-arm 
selection. The model-dependence of the acceptance corrections results in an uncertainty for 
eacc(iTiodels) of 0.5% (1.6%) for the two selections, respectively. 

Since less than half of the diffractive events, mostly with a high-mass diffractive system, pass 
the hadron-level selection, the uncertainty on Cvis is smaller than that on Cacc- The 1 standard 
deviation differences found among the four hadronic interaction models on the hadron-level 
correction, Cvis/ propagate into uncertainties on cTVis of 0.4% (0.2%) for the single-arm (double¬ 
arm) selection. The subtraction of photon-proton events (with the N.yp uncertainty shown in 
Table [^, results in an uncertainty of 0.6% (<0.1%) on (T-nei and PVis. The uncertainty on Nnoise 
propagates into a 1.3% (0.2%) uncertainty in the final cross sections. The effect on the event 
selection of the radiation damage in the HF fibres is assessed by rescaling the signals of the 
simulated HF response to match data in segments of pseudorapidity. The rescaling factors are 
calculated using the average response produced by EPOS, HIJING, and QGSJETIL These scaling 
factors are found to be consistent with the observed radiation damage of HF and range from 
1 to 0.67, depending on pseudorapidity. The amount of radiation damage is estimated from 
a comparison of dE/df] distributions measured in proton-proton collisions at \/s = 2.76 TeV 
recorded in 2013 and in 2010. The systematic uncertainty induced on the cross section by this 
approach is estimated by repeating the measurement without the radiation damage correction, 
which introduces an effect of 1.7% (0.8%) on the cross section. As a further check of the HF 
tower energy resolution, the cross sections are computed by increasing the selection thresholds 
to Ehf > lOGeV (5GeV). To account for both effects, a systematic uncertainty on the cross 
section of 0.6% (0.4%) is added. The cross sections measured for the two beam directions are 
found to be consistent. Gonsequently, no dedicated systematic uncertainty is assigned to this 
effect. 

All the different sources of uncertainty of the measurement are listed in Table for the single¬ 
arm and double-arm event selections. The three derived cross sections have different system¬ 
atic uncertainties since not all contributions are relevant to each of them. For Pinei/ all un¬ 
certainties but the one due to the hadron-level correction contribute. The total uncorrelated 
systematic uncertainty is therefore 2.5% (2.2%) for the single-arm (double-arm) selection. For 
(Tyis, the dominant uncertainty is due to the hadron-level correction instead of the correction 
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Table 2: List of the systematic uncertainties, propagated into the final pPb cross sections, for 
the two event selections. 


Source of uncertainty 

Single-arm 

Double-arm 

Noise subtraction (Nnoise) 

1.3% 

0.2% 

Pileup correction (/pu) 

<0.1% 

<0.1% 

Acceptance (Cacc(models)) 

0.5% 

1.6% 

Acceptance (eacc(Giiff)) 

0.8% 

1.1% 

Hadron-level correction (Cvis) 

0.4% 

0.2% 

Photon-proton subtraction (N^p) 

0.6% 

<0.1% 

Detector simulation 

1.7% 

0.8% 

HP energy thresholds 

0.6% 

0.4% 

Integrated luminosity (C) 

3.5% 

3.5% 


for Cacc- The value of the uncertainty is therefore reduced to 2.3% (0.9%). The uncertainties 
for detector simulation and photon-proton correction do not contribute to cTo^g and, hence, its 
uncertainty becomes 1.4% (0.5%). For all cross sections, a (dominant) integrated luminosity un¬ 
certainty of 3.5% is added. The main contributions to the latter arise from the model used to 
describe the beam profile, the length scale of the beam displacement, and the bunch-to-bunch 
variations 1321 . 


4 Results and summary 

The measured cross sections for the single-arm and double-arm event selections are listed in 
Table compared to the predictions of the hadronic interaction models dpmjet, EPOS, and 
QGSJETll. Due to the different acceptance, the extrapolations from the hadron-level to the in¬ 
elastic cross section are of different magnitude, but the models reproduce well the approxi¬ 
mately 65 mb difference between the two selections. The values of the inelastic cross sections 
obtained from the single-arm and double-arm methods agree well within the uncertainties. 

Table 3: Summary of cross sections obtained from the two different event selections. The accep¬ 
tance definition for PVig is based on the production of stable particles within 3 < \rj\ < 5 with 
momentum phf > 21.3GeV/c (11.3GeV/c) for the single-arm (double-arm) event selections. 



Selection 

G)bs (mb) 

Gvis (mb) 

Gnel (mb) 

Data 

Single-arm 

2003±76 

1937±82 

2063±89 

Double-arm 

1873±66 

1872 ±68 

2059±85 

EPOS-LHC 

Single-arm 

Double-arm 

— 

1947 

1883 

2082 

QGSJETII-04 

Single-arm 

Double-arm 

— 

2059 

1998 

2181 

DPMJET 3.06 

Single-arm 

Double-arm 

— 

2116 

2055 

2166 


The final Uinei value is obtained by taking the weighted average of the measured values in the 
two event selections. The statistical uncertainties and the uncertainty on the luminosity are 
correlated between the selections. The degree of correlation among the remaining systematic 
uncertainties is much smaller and they are taken as uncorrelated. This yields a final result for 
the inelastic hadronic cross section of 


Gnei(pTT) = 2061 ± 3 (stat) ± 34 (syst) ± 72 (lumi) mb. 
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4 Results and summary 



Figure 3: Inelastic hadronic cross sections for pPb collisions as a function of the centre-of-mass 
energy. The measurement described here (circle, with error bars obtained from the quadratic 
sum of all uncertainties) is compared to lower energy data (squares and triangles) B |39l 
and to different model predictions (curves). 

This result is shown in Fig. [^compared to other measurements at different centre-of-mass en¬ 
ergies and to various theoretical predictions. A pPb cross section was also measured by the 
ALICE Collaboration, amounting to 2090-2120 mb with an uncertainty of 70 mb |4TJ. A di¬ 
rect comparison of this observed cross section to the one measured in the present analysis is 
not possible due to the unknown to us ALICE detector acceptance and possible contamination 
from noise and photon-proton interactions. 

The inelastic cross section measured by the CMS experiment is compared to the Glauber-model 
prediction (solid curve in Fig.|^ obtained using a pp inelastic cross section at \/s = 5.02 TeV of 
70.0 ±1.5 mb, derived from the COMPETE parametrisation 1421 including the measurement of 
the TOTEM Collaboration at = 7 TeV Il43l (where the assigned uncertainty is that measured 
by the latter). The Glauber calculation yields 2130 ± 40 mb and is compatible with the mea¬ 
surement presented here indicating that effects neglected by the calculation (such as nucleon 
correlations and screening) are either small or approximately cancel out. The experimental re¬ 
sult is also consistent with the prediction of the DIPSY model Il44l l45l based on a dipole-model 
approach including parton saturation and multiple-scattering. Among the Gribov-Regge mod¬ 
els, the EPOS prediction is compatible with the measurement within uncertainties, whereas DP- 
MJET and QGSJETll predict a value more than 1 standard deviation above the data, with a larger 
discrepancy appearing for the Gvis cross sections (Table |^. The EPOS and QGSJETll models are 
commonly used for cosmic ray air shower simulations. Thus, at the corresponding cosmic ray 
proton energies of Ea = s/ (2mp) = eV, where nip is the mass of the proton, there are 
no indications for data-model deviations above «5% in the proton-lead collisions studied here. 
Note, however, that our measurement deals with an ion much heavier than those involved 
in proton-air interactions. Gorrections to the Glauber model are possibly larger in the latter 
case 13 SSI. In summary, the measurement of the cross sections in pPb collisions presented 
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here is the first such fully corrected measurement at multi-TeV energies and, thus, provides 
important constraints on hadronic interaction models commonly used in high-energy heavy 
ion and cosmic ray physics. 
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